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ABSTRACT 

The dynamics of globular cluster systems (GCSs) around galaxies are often used to 
assess the total enclosed mass, and even to constrain the dark matter distribution. The 
globular cluster system of a galaxy is typically assumed to be in dynamical equilibrium 
within the potential of the host galaxy. However cluster galaxies are subjected to a 
rapidly evolving and, at times, violently destructive tidal field. We investigate the 
impact of the harassment on the dynamics of GCs surrounding early type cluster 
dwarfs, using numerical simulations. We find that the dynamical behaviour of the 
GCS is strongly influenced by the fraction of bound dark matter /dm remaining in 
the galaxy. Only when /dm falls to ~ 15%, do stars and GCs begin to be stripped. Still 
the observed GC velocity dispersion can be used to measure the true enclosed mass to 
within a factor of 2, even when /dm falls as low as ~ 3%. This is possible partly because 
unbound GCs quickly separate from the galaxy body. However even the distribution 
of bound GCs may spatially expand by a factor of 2-3. Once /dm falls into the < 3% 
regime, the galaxy is close to complete disruption, and GCS dynamics can no longer 
be used to reliably estimate the enclosed mass. In this regime, the remaining bound 
GCS may spatially expand by a factor of 4 to 8. It may be possible to test if a galaxy 
is in this regime by measuring the dynamics of the stellar disk. We demonstrate that 
if a stellar disk is rotationally supported, it is likely that a galaxy has sufficient dark 
matter, that the dynamics of the GCS can be used to reliably estimate the enclosed 
mass. 

Key words: methods: N-body simulations — galaxies: clusters: general — galaxies: 
evolution — galaxies: kinematics and dynamics — galaxies: star clusters: general — 
galaxies: dwarf 



1 INTRODUCTION 

1.1 Cluster galaxies in the hierarchical galaxy 
formation scenario 

In the current framework of structure formation, cold dark 
matter (CDM) haloes grow hierarchically by the merging of 
smaller, virialised haloes. After hifall, these subhaloes con- 
tinue orbiting within the main halo, with tidal interactions 
stripping their masses from the outside-in and, if massive 
enough, spiraling in as dynamical friction slows down their 
velocities. Thus, depending on their concentration and or- 
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bital distribution, subhaloes can lose substantial mass and 
eventually disrupt or merge at the halo centre. This is a 
process which is specially efficient in the inner, dense cen- 
tral regions of the main halo. The interplay between contin- 
uous accretion and disruption results in haloes containing 
~ 10-20% of their mass in the form of self-bound subhaloes 
jGhigna et al.lll998l : iGill et alj|2004 ). Galaxies form in DM 
haloes as baryons cool and contract towards the centre of 
the halo potential well (|White k, Reed [T978t ) . In this con- 
text, as smaller haloes are accreted and become subhaloes 
of larger hosts, their galaxies can be identified as a satellite 
population. Virialised galaxy clusters conform the paradigm 
of this scenario, with several hundreds of galaxies orbiting 
within 10 14 — 10 15 Mq haloes. Close, high speed tidal encoun- 
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ters between cluster galaxies are predicted to be frequent, 
and can result i n substantial mass loss - a process known 
as 'harassment' (jMoore et al.lll99rjh . Given that galaxies re- 
side deep within their halo potential well, the efficiencies 
of D M and baryons strippin g are not straightforwardly re- 
lated (jPefiarrubia et aD2008l ). and understanding the fate of 
satellite galaxies has become on e of the greatest challe nges 
for models of galaxy formation (jWetzel fc Whitell2010l ). 

The current situation is somewhat paradoxical. While 
theory and numerical simulations provide accurate descrip- 
tions of the stripping process and mass evolution of DM 
subhaloes, the evolution of their baryonic component is, 
in comparison, more difficult to track. For instance, most 
semianalytic models have, so far, simply assumed that the 
stellar mass remains intact until satellites eventually merge 
with their central galaxy |De Lucia fc BlaizotJ l2007h . The 
opposite is true from an observational point of view: while 
the stellar content of cluster galaxies is routinely measured, 
very little is known about their true DM content and the 
properties of their haloes. The reason for this discrepancy is 
that traditional mass determination methods face significant 
challenges when it comes to determining the total DM halo 
mass. Stellar dynamics studies are very expensive in terms 
of telescope time, and therefore have been often limited to 
the inner (< r e ) regions of luminous galaxies - which are al- 
ready known to be not DM-dominated. Only recently have 
observations been able to probe regions as far as ~ 3r e , but 
the stellar component sti ll provides the major contribution 
to th e system dynamics (|Proctor et al.l [20091 ; ISpolaor et alj 
I2010T I. The situation is even more complex in the case of 
low-mass cluster galaxies, where the intrinsically low surface 
brightness restrict stellar kinematics measurements to the 
inner effective radius (e.g.. [Chilingarianll2009l ; iToloba et alj 
l201ll ). Gravitational lensing has been suggested to be the 
most reliable method for determining the extent and total 
mass of cluster subhaloes, but the induced shear is so weak 
that the approach br eaks down at the exp ected masses of 
dwarf- hosting haloes l|Limousin et al.ll2009l ). 



1.2 Globular cluster systems as dynamical probes 

Globular cluster systems (GCSs) perhaps provide the only 
means by which the kinematics of cluster galaxies at large 
galactocentric radii can be probed. Globular clusters are 
old, relatively luminous systems exhibiting very extended 
galactocentric distributions. Furthermore, and contrary to 
planetary nebulae - the other stellar population tradition- 
ally used as a kin ematic tracer of the host halo (e.g., 
ICoccato et alll2009] ) - GCs are found in abundance in clus- 
ter galaxies across a wide stellar mass range. They are 
rather populous systems in all Brightest Cluster Galax- 
ies, but are found as well in cluster early-type dwarfs 
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Georgiev et al.ll2009l ; [Georgiev et aTll2010l ). They are known 
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to be effective kinematic probes for luminous galaxies, 
where GCSs have been used to obtain rotation curves 
and radial velocity dispersion profiles extendi ng well be- 
yond r > 3 r e for about a do zen galaxies JPuzia et alj 
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the consistency between these results and those derived 
from integrated stellar kinematics in the overlapping regions 
strongly argues in favour of GCSs being true tracers of the 
old, metal-poor stellar haloes in massive galaxies. 

The low-mass end of the stellar mass function is a 
much lesser explored territory. Perhaps surprisingly, it is now 
well known that the family of luminous Vir go dEs (M* > 
2 x 10 8 M B ) tend to co ntain populous GCSs (|MilIer fc Lotzl 
120071 ; IPeng et al.ll2008l ). The majority of these massive, GC- 
rich dEs are found within the nucleated subclass, and their 
high GC mass specific frequencies represent a challenge: hav- 
ing spent a considerable amount of time within the harsh 
cluster potential, it is intriguing how have these systems 
managed to retain such a high abundance of GCs. The first 
kinematical study of GCs in a Virgo cluster dwarf galaxy 
was that of the dE VCC 1087 bv lBeaslev et al.l (120061 ). Spec- 
troscopy of 12 GCs out to r « 6 r e in this « 3 x 10 9 M o 
dE galaxy surprisingly show evidence of dynamically sig- 
nificant rotation in the GCS, with (v/a)* > 1.6. The im- 
plied dynamical mass is rel atively high (Mdyn/M j, ~ 8), 
a result later confirmed by iBeaslev et al.l (|2009l ) for two 
other luminous Virgo dEsQ Interestingly, the bulk of the 
dEs stellar bodies do not seem to share the GCS rotation, 
and these authors suggest that these features might be con- 
sistent with evolution through harassment from a late-type 
disc g alaxy. On the other hand, ISanchez-Janssen fc Aguerril 
((2OI2} ) point out that the GC abundance of these dEs is 
too high (factors 2-8) compared to the one expected from 
such progenitors, and argue that high total masses might 
be required for dEs to maintain these populous GCSs dur- 
ing their orbital evolution within Virgo. It thus appears 
that kinematical studies of GCSs in dwarf galaxies offer a 
unique opportunity to determine their total mass content 
at large radii - a region currently inaccessible to any other 
technique and where DM is expected to be dominant, in 
agreement with those findings. They can help us bridge the 
gap between the very high Mdyn/Mj, rati os characteristic 
of Local Group dSphs (|Strigari et al.l [2*0081 ) and the rather 
low values (Md yn /M± ~ 3) systematical ly found from stel- 
lar kinematics of c luster and group dEs (|Toloba et al.ll201ll ; 
iForbes et a"i1l201ll ). 

However, as previously stated, the orbital evolution 
of subhaloes within clusters can significantly influence the 
dynamical status of GCSs. For instance, cluster galaxies 
are subjected to a rapidly evolving and, at times, vio- 
lently destructive tidal field that can bring the system 
out of dynamical equilibrium - an implicit assumption in 
all dynamical models. Also high-speed tidal encounters in 
a cluster are expected to result in significant dark mat- 



ter mass loss in dwarf galaxies dMastropietro et al.l 120051 : 
lAguerri fc Gonzalez-Garci'al 120091 : ISmith et all l2oiof ). and 



20091 : ISchuberth et al.ll201fj| : iFoster et ; all 12011! ). Moreover" 



the GCS must settle within the new mass distribution. More- 
over, the tidal evolution within the cluster can result in some 
GCs being stripped, no longer tracing the original host halo 
and biasing results if included within observational samples. 
All these effects may significantly impact the implied 



1 Following previous work (e.g., Boselli et al. 2008), with this 
definition we refer to any low-mass, gas-free and non-star-forming 
cluster galaxy. Throughout this paper they will all be termed dE 
for simplicity. 
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dE dynamical masses and yet, ironically, harassment is ex- 
pected to result in substantial dark matter mass loss that 
may well be detectable by observing GCS dynamics. It is 
therefore fundamental that the impact of fast close encoun- 
ters on mass estimates from GCSs is well understood. 

Here we provide such an investigation by studying the 
impact of cluster harassment on the dynamics of GCSs sur- 
rounding dE galaxies, using numerical simulations. In this 
paper, we focus on GCS dynamics as a means of measur- 
ing the enclosed mass. In a future publication we shall fully 
investigate additional topics such as induced rotation, spa- 
tial distribution, and the fate of the stripped GCs. We de- 
scribe our numerical code, the cluster harassment model, 
and galaxy model setup in Section 2, we present our results 
for mass-loss, effects on spatial distribution, and dynamics 
of the GCSs in Section 3, and we summarise and draw con- 
clusions in Section 4. 



2 SETUP 



2.1 The code 



In this study we m ake use of 'gf' (|Williams fc Nelson! 
l200ll : IWilliams1ll998h . which is a Treecode-SPH algorithm 
that operates primarily using the techniques described in 
iHernauist fc Kata (jT989[) - The dE galaxies we consider are 
gas-free and non-star-forming. We therefore do not include 
an SPH component or star formation recipe to our models, 
and the code operates purely as a gravitational code with- 
out considering hydrodynamics, 'gf' has been parallelised 
to operate simultaneously on multiple processors to decrease 
simulation run-times. The Treecode allows for rapid calcu- 
lation of gravitational accelerations. In all simulations, the 
gravitational softening length, e, is fixed for all particles at 
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in common with the harassment simu- 



lations of iMastropietro et al.l (|2005h . Gravitational acceler- 
ations are evaluated to quadrupole order, using an open- 
ing angle 8 C — 0.7. A second order individual particle 
timestep scheme w as utilised to improve effic iency following 
the methodology of lHernquist fc Katd (1 19891 ). Each particle 
was assigned a time-step that is a power of two division of 
the simulation block timestep, with a minimum timestep of 
~0.5 yrs. Assignment of time-st eps fo r collis ionless particles 
is controlled by the criteri a of Katj ( 199ll) . For details of 
code testing, please refer to I Williams! ( 19981 ). 



2.2 Harassment model 

We model the tidal field of a Virgo-like cluster using a two 
component model; the first is a dynamical tidal field associ- 
ated with individual 'harasser' galaxies, and the second is a 
static tidal field to represent the rest of the mass of the clus- 
ter that is not associated with galaxies - referred to as the 
background cluster potential herein. The potential field of 
each individual harasser and the background is represented 
using the analytical for m for an NFW density distribution 
(|Lokas fc Mamonll200j ): 



where g c — l/[ln(l + c) — c/(l + c)], c is the concentration 
parameter, r s is a characteristic radial scalelength, and M200 
is the virial mass. 

We fully define the background cluster potential by 
choosing its virial mass M2oo=l-6x 1O 14 M0, concentra- 
tion c=4 and virial radius r2oo=H00 kpc. This choice of 
m ass is in close agreem ent with the Virgo cluster model 
of IVollmer et al.l (|200ll ) , and the concentration is typical 
for c luster-mass halos in CDM simulations (e.g. I Gill et all 
l2004f ). Now we give each harasser galaxy its own individual 
NFW potential field. To make this field dynamical, a one-off 
N-body simulation of a cluster-mass halo is conducted, and 
the time-evolving coordinates of a fraction of the particles 
are logged. This log now provides the position-evolution of 
the centre of each of the harassing haloes, where each indi- 
vidual galaxy potential field will be super-imposed. We re- 
quire a mass and concentration for each harasser's potential 
field to be defin ed. For mass, we use the Sch echter function 
Schechterl 19761 ) with parameters provided in lSandage et all 



1989) , and assume a constant mass-to-light ratio: 
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where a = -1.25, a nd (M/L) K = 20. T he fitting pa- 
rameters provided in ISandage et alj (Il985l) are based on 
the Virgo cluster catalogue (|Binggeli et"afl Il985l ). and are 
therefore complete down to a minimum galaxy luminosity 
of ~ 2.5 x 10 7 L(7) assum ing a distance modulus of 31.0 
for Virgo jMei et al.ll2007l ). We resolve all harasser galax- 
ies down to one-hundredth the mass of our standard dwarf 
model (1O 9 M0). Completing the integral in Equation [2] pro- 
duces a total mass in galaxies that is 14.5% of the total clus- 
ter mass. This is in reasonable agreement with the ~ 10% 
percentage of a cluster's mass lock ed up in sub-str ucture, 
found in ACDM simulations - see iGill et ail (|2004l ). This 
process produces 733 harassing galaxies in total. The vast 
majority of harasser galaxies are dwarfs, with ~ 90% of 
the total number of harassers having a mass less than the 
standard dwarf galaxy model. The maximum mass galax- 
ies produced are ~ 1O 12 M0. The value of the mass-to- 
light ratio chosen is consistent, if not an upper lim it for 
galaxies in this mass-range (see iGilmore et al.l 120071 ) . The 
concentration c of each halo is chosen to follow the trend 
of higher concentration for lower mass objects as found in 
ACDM simulations. This is a chieved by a fi t to co ncentra- 
tion values with mass found in lNavarro et al.l (|l996l ) produc- 
ing c = -31og(M ha i o (M )) + 52. This fit produces c ~ 15 
for a 6 x 10 11 M mass halo, and c ~ 20 for a 1 x 10 10 
Mq mass halo. The total mass of our model cluster is fixed, 
whereas real clusters tend to increase their mass by accre- 
tion. Therefore we study the effects of harassment over a 
limited duration of 2.5 Gyr. In this way, the inaccuracy in- 
troduced by assuming a fixed cluster mass is also limited. 
Further properties of the harassment model can be found in 
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ISmith et al1(|2010l ). 



* = -g c GM 20 o 



ln(l + (r/r s )) 



(1) 



2.3 Galaxy models 

Our disk galaxy models consist of 3 components ; an NFW 
dark matter halo (|Navarro. Frenk. fc Whitelll996l ). an expo- 
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nential disk of stars, and a spheroidal distribution of globular 
cluster particles following a Hernquist profile. 

2.3.1 The dark matter halo 

The dark matter halo of the disk galaxy model has an NFW 
density profile. The NFW profile has the form: 



p(r) 



(~)(i + £) 2 



(3) 



where po is the central density. Given that the NFW model 
has a divergent total mass, the profile is typically truncated 
outside the Virial radius, r2oo = r s c. 

Positions and velocities are assigned to the dark mat- 
ter particles using the pub lically available algorithm mkhalo 
from the nemo repository (|McMillan fe Dehnenll2007r ). Dark 
matter halos produced in this manner are evolved in isola- 
tion for 2.5 Gyr to test stability, and are found to be highly 
stable. 

Our standard model has a dark matter halo mass of 
fO n Mo, consisting of f 00, 000 dark matter particles, with a 
concentration c = f 4. The Virial radius is r2oo = 95 kpc. The 
peak circular velocity of the halo is 88 km s _1 at a radius of 
15 kpc. If the halo were placed at a radius of 200 kpc within 
the smooth background potential of the cluster, the ratio 
of its central density is ~ 5 times that of the surrounding 
ambient medium. 



2.3.2 The disk 

The stellar disk of the galaxy has an exponential form 



£(7?) = £ exp(i?/i? d ) 



(4) 



where E is the surface density, Eo is central surface density, 
R is radius within the disk, and Rd is the scalelength of the 
disk. 

The scalelength and mass of the stellar disk is chosen 
to approximately match the observed properties of luminous 
Virgo dEs. These are characterised by close-to-exponential 
luminosity profiles w ith effective radii of the order of 1 kpc 
|janz fc Lisker|[2008T l. Our fiducial model has a total stella r 
mass of 3.0 x 10 9 M Q (3% of the halo mass; |Peng et al.ll2008h . 
and is formed from 20,000 star particles, distributed with 
i?d = 0.75 kpc, resulting in an effective radius of ~1.25 kpc. 

A radially varying velocity d ispersion is ch osen that en- 
sures the disk is Toomre stable |Toomr eHl964T l at all radii. 
In practice, a Toomre parameter of Q > 1.5 is required 
throughout the stellar disk to ensure stability. We create 
a dispersion dominated disk with v ro t/vdi sp ~ 0.25. This 
results in a thick disk with axial ratio b/a ~ 0.6 in agree- 
ment with t he observed axial ratios of dwarf galaxies of thi s 
luminosity |Lisker et alj|2007l ; ISanchez- Janssen et al.ll2010t ). 
A full description of the pr ocedures followed t o set up the 
galaxy disk can be found in ISmith et all |2010T l. 

2.3.3 The globular cluster distribution 

Each globular cluster is treated as a single N-body parti- 
cle. The standard model has 60 equal-mass globular cluster 
particles with a summed total mass of 1.8 x 10 7 M , corre- 
sponding to a specific stellar mass fraction of 0.6 per cent 
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Figure 1. Surface density profiles of the dark matter (solid line), 
stars (dashed line), and GCs (dot-dashed line), for all the galaxy 
models. 



(|Peng et al.ll2008l : lGeorgiev et al.ll20ld ). We model the glob- 
ular cluster particle distribution as a Hernquist sphere: 



p(r) 



M h r h 



2n r (r + rh) 3 



(5) 



where Mh is the total mass of the Hernquist sphere, rh is 
the Hernquist scalelength, and r is radius. 

The distribution is radially truncated at a cut-off ra- 
dius of 7.5 kpc. The observed GCS density distribution in 
luminous dEs depends on projected ra dius roughly as R~ 
jBeaslev et al.l 120061 : IPuzia et al.ll2004l ). We obtain an ad- 
equate match to this radial dependency by choosing the 
Hernquist scalelength rh = 3.75 kpc (half the cut-off radius). 
Particle velocities are assigned using the Jeans equation for 
an isotropic dispersion supported system. 

Finally, we note that all models (combined dark matter 
halo, stellar disk and surrounding GCs) are evolved in isola- 
tion for 2.5 Gyr to ensure stability, before introduction into 
the cluster environment. 
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2.4 A parameter study 

Our standard model is meant to roughly approximate the 
subclass of luminous dEs - the only ones that have had 
their GCSs studied in sufficient detail. However, we are in- 
terested in understanding how the effects of harassment on 
any galaxy are dependent on the parameters of that partic- 
ular galaxy model. We systematically vary one parameter at 
a time creating a set of 8 galaxy models as listed below. 
Halo concentration: We test a high concentration halo model 
('High c') with c=30, and a low concentration halo model 
('Low c') with c=5. The standard model has c=14. 
Halo mass: We test a high mass halo model ('High Mass') 
with M200 = 3.0 x 10 Mq, and a low mass halo model ('Low 
Mass') with M200 = 0.333 x 1O 11 M . The standard model 
has M200 = 1.0 x lO n M . Measured within 7.5 kpc (6 ef- 
fective radii) , the dynamical mass-to-light ratio for the high 
mass, standard, and low mass halo model is Mdyn/M* ^2 9, 
5, and 4 respectively. 

Globular cluster distribution: We test a concentrated GCS 
('Cone GCS') with Hernquist scalelength rh = 1.875 kpc, 
truncated at 3.75 kpc, and an extended GCS ('Ext GCS') 
with rh = 7.5 kpc, truncated at 15 kpc. For comparison, 
the standard model has rh = 3.75 kpc, and is truncated at 

7.5 kpc. 

Disk scalelength: We test a small stellar disk model ('Small 
Disk') with exponential scalelength rd = 0.375 kpc, whereas 
the standard model has r^ = 0.75 kpc. The surface density 
profile of each model is shown in Figure [T] 

We note that we do not attempt to simulate fainter 
cluster dEs or dSphs. The strong dependence of GC specific 
frequency on stellar m ass for these galaxies (|Miller fc Lot j 
l2007l ; |p"e"ng et al-lfeoOct ) imply that essentially all dwarfs hav- 
ing M t <3x 1O 8 M host less than 10 metal-poor GCs. We 
therefore expect that the range of parameter space we cover 
should encompass the properties of most early type dwarfs 
who contain sufficient numbers of GCs to be useful for an 
enclosed mass measurement. 



2.5 Orbits 

Our harassment model enables us to evolve a chosen galaxy 
model along a fixed orbit through the cluster. Each galaxy 
model then experiences an identical sequence of tidal en- 
counters. Thereby we can fairly compare the responses of 
different galaxy models to identical harassment histories. If 
the galaxy takes an alternative trajectory through the clus- 
ter, harassment histories are highly stochastic resulting in 
widely differing numbers and strengths of tidal encounters 
along the orbit. 

Although stochastic from orbit to orbit, stronger ha- 
rassment typically occurs for orbits that spend more time 
near the cluster centre where the number-density of harasser 
galaxies is highest, and the tid al force from the smooth 
background potential is largest jSmith et alj|20icl ). In or- 
der to sample a range of harassment histories from strong 
to weak, we conduct 16 different orbits. The parameters of 
these orbits are shown in Table[T] The first eight orbits ('or- 
bit P1-P8') are on a plunging orbit with an initial radius 
of ~ 987 kpc. The initial velocity is 450 km s _1 and purely 
tangential. In the absence of harassers, orbit P1-P8 would 
have a pericentre of ~ 230 kpc. The latter eight orbits ('or- 



Orbit 


Orbit 


I'init 


I'm in 


Tmax 


1 ' peri 


foM 


f* fee 


number 


type 


(kpc) (kpc) (kpc) 




(standard) 


PI 


Plunge 




987 


y / 


1 




8Q Q8 


P2 


Plunge 


987 


201 


1250 


1 


0.02 


U.'-tO u.oo 


P3 


Plunge 


Jo / 


ZOO 


1024 


1 


u.zo 


n Q7 1 on 


PI 


Plunge 


987 


235 


1793 


1 


0.18 


1.00 0.98 


P5 


Plunge 


987 


309 


973 


1 


0.23 


0.98 1.00 


P6 


Plunge 


987 


313 


973 


1 


0.35 


0.98 1.00 


P7 


Plunge 


987 


318 


972 


1 


0.01 


0.22 0.20 


PS 


Plunge 


987 


200 


1953 


1 


0.13 


1.00 0.97 


CI 


Circular 


200 


102 


332 


3 


0.00 


0.20 0.20 


02 


Circular 


200 


111 


101 


3 


0.04 


0.68 0.70 


C3 


Circular 


200 


110 


197 


3 


0.02 


0.57 0.47 


01 


Circular 


200 


124 


201 


3 


0.00 


0.00 0.00 


05 


Circular 


200 


52 


334 


3 


0.00 


0.00 0.00 


06 


Circular 


200 


124 


240 


3 


0.08 


0.86 0.92 


07 


Circular 


200 


92 


232 


1 


0.00 


0.00 0.00 


08 


Circular 


200 


107 


292 


3 


0.01 


0.30 0.23 



Table 1. Orbit parameters. Columns (from left to right): (1) 
Orbit label, (2) Orbit type (plunging or circular), (3) Initial radius 
riniti (4) Minimum radius in cluster r m ; n , (5) Maximum radius 
in cluster r max , (6) Number of pericentre passes N pGr i over 2.5 
Gyr, bound mass fraction of the standard model after 2.5 Gyr 
of harassment for the (7) dark matter foMi (8) stars f*, and (9) 
globular clusters {qq 

bit C1-C8') have an initial radius of only 200 kpc from the 
cluster centre. The initial velocity is 1500 km s _1 and purely 
tangential. In the absence of harassers, orbits C1-C8 would 
be circular. However, harassment causes all the orbits to de- 
viate, resulting in variation in the minimum and maximum 
radius from the cluster centre (see columns labelled 'r m i n ' 
and 'r max '). 

The column labelled 'N pcr i' indicates the number of 
times the galaxy passes pericentre along its orbit. The fi- 
nal three columns (labelled /dm, /*, and /gc) are the final 
bound fraction of the dark matter, stars and globular clus- 
ters, respectively, for the standard model. N pcr i, /dm, /* and 
/gc are measured after 2.5 Gyr of harassment along that or- 
bit, /dm can be considered an indication of the strength 
of the harassment along that orbit (i.e. a low /dm indi- 
cates large mass loss of dark matter through strong harass- 
ment). Orbit P1-P8 typically experience weak harassment, 
as they make only a single plunging pass of their pericen- 
tre, and so spend little time near the cluster centre where 
the highest number density of harassers can be found. How- 
ever strong mass loss can still occur in chance strong tidal 
encounters (e.g. see orbit P2). Orbits C1-C8 typically expe- 
rience strong harassment as they remain close to the cluster 
centre throughout their orbit. In total, we conduct a total 
of 16 orbits on 8 galaxy models, resulting in a total of 128 
harassment simulations. 



2.6 Measuring bound fractions 

The following technique is found to be a reliable, and repro- 
ducable method for measuring the fraction of particles that 
remain bound to the galaxy at any instant. In the first step, 
all particles within a 2.5 kpc radius of the centre of den- 
sity of the galaxy are selected. In practice, the final bound 
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Figure 2. Plots of the time evolution of the fraction of bound 
particles of dark matter (black solid line), stars (red dashed line), 
and GCs (green dash-dotted line) for the standard galaxy model, 
undergoing harassment. The left panel is orbit P2 and is a medium 
strength harassment case with roughly half the stars and GCs 
stripped. The right panel is orbit P8 which is a weak harassment 
case, with ~ 10% of stars and GCs lost 

fractions are found to be very insensitive to this choice of 
radius. Each particle confirms if it is bound to the other par- 
ticles within this radius. Those that are bound to each other 
are considered a bound core within the galaxy center. Now 
an iterative procedure begins. In each iteration, all particles 
are tested to see if they are bound to the bound core. If 
they are, then their mass is added to the bound core. We 
call this method of growing the mass of the bound core the 
'snowballing' method. Then the iteration is repeated until 
the bound mass of the galaxy increases by no more than 1% 
between iterations. In practice the total bound mass is found 
within 5-10 iterations. We find the snowballing method of 
measuring bound masses to be robust and trustworthy. A 
full description and test ing of this approa ch will be pre- 
sented in a future paper (jFarias et al.ll2012l ; in prep.) 



3 RESULTS 

3.1 Representative examples 

In Figure [2] we show two examples of mass loss of dark mat- 
ter, stars and GCs, as a result of harassment acting on the 
standard model galaxy. The left hand panel is orbit P2 - 
an intermediate strength harassment case. The right hand 
panel is orbit P8, which is a case of weak harassment. How- 
ever in both substantial fraction (> 80%) of the 
dark matter is unbound. This occurs rather gradually un- 
der the influence of the smooth, background potential of the 
cluster (from mass not associated with individual galaxies). 
However it can also occur rather rapidly when a high speed 
encounter with another cluster galaxy occurs (e.g. see the 
left panel at t=1.4 Gyr). In both the left and right panel, 
we note that the stars and globular clusters do not begin to 
be stripped until the dark matter fraction is low (~10-20%). 
We shall better quantify this value in section 13.41 

In all our simulations, we treat each individual GC as a 
point mass, whereas in reality an ind ividual GC can be sub- 
ject to tidally induced mass loss (e.g. iBaumgardt fc Kroupal 
120071 : ICieles et al1l201ll : iKiipper et al.ll2012f l . Assuming that 
each of our GC contains 10 s M© within a radius of 3 pc 
(equivalent to NGC104 or similar), we can calculate the 
Roche limit of a GC. The Roche limit is the maximum tidal 



force the GC can experience without being disrupted. We 
then monitor the tidal forces that each GC experiences dur- 
ing the orbit CI infall. Orbit CI is a case of strong harass- 
ment. In this simulation strong harassment results in heavy 
stripping of dark matter (/dm = 0.00), stars (/* = 0.20), 
and GCs (/gc = 0.20). However, we find that GCs do not 
experience tidal forces greater than 0.1% of their own Roche 
limit, over the total 2.5 Gyr of harassment. Although many 
are stripped from the their host galaxy, they would not be 
completely disrupted by the tidal forces. Mass loss from real 
GCs also occurs due to internal processes operating within 
the G C, for examp le two-body encounters and stellar evo- 
lution (|Gielesll2012r i. However, it is unlikely that these pro- 
cesses would result in disruption of many GCs during the 
short duration of the harassment simulations (2.5 Gyr). As 
our simulations use softened gravity between particles, we 
suppress the effects of dynamical friction acting on our GCs, 
whi ch could result in their sinking towards the galaxy cen- 
tre (|Oh fc LinllioOtil) . However we expect the effect to be 
negligible for our model galaxy (i.e. our higher mass galaxy, 
with lower mass GCs results in a dynamical friction drag 
force that is ^2 -3 orders of magnitudes weaker than in the 
lOh fc LinllioOOl study). 

3.2 Mass- loss dependency on galaxy properties 

We quantify the harassment induced mass loss that occurs 
to the dark matter, stars and GCs of all 8 galaxy models, for 
the 16 orbits we consider. To do this, we measure the fraction 
of dark matter (/dm), stars (/*) and globular clusters (/gc) 
that remain bound to the galaxy, after 2.5 Gyr in the clus- 
ter environment. Although we only present the bound mass 
fractions at 2.5 Gyr, we find the relative mass loss between 
models at earlier times (e.g. at 2 Gyr) is fairly represented. 
A low bound fraction therefore reflects that lots of mass has 
been stripped. We additionally conduct simulations where 
all the harasser galaxies are removed, and mass loss occurs 
purely as a result of the smooth backrgound potential of the 
cluster. We refer to these simulations as 'no-harasser' runs, 
and conduct them for all model galaxies, on the plunging 
and circular orbits. The results are shown in Figure [3] 

In the first, second, and third row we plot /dm, /*, 
and /gc respectively, against the orbit type (plunging or 
circular). In the bottom panel the ratio of /gc//* is plot- 
ted - this ratio can indicate if stars or globular clusters 
have been preferentially lost, and is directly proportional 
to the GC mass s pecifi c frequency, T N = jV S c/(M*/10 9 M Q ) 
()Zepf fc Ashman1ll993h . Columns indicate the model galaxy, 
as indicated by the column header. Plunging orbits are indi- 
cated by symbols with orange-to-red colours. Circular orbits 
are indicated by symbols with light-to-dark blue colours. A 
particular orbit (e.g. orbit CI) has the same colour symbol 
in every panel. Circular symbols are results from the ha- 
rassment simulations. A star symbol indicates a no-harasser 
orbit, where mass loss is purely from the tidal forces of the 
smooth background potential of the cluster. 

First, comparing the upper row to second row, it is clear 
that dark matter is always preferentially stripped in com- 
parison to stars. The dark matter is more extended and so 
more easily unbound. However if we had instead plotted the 
number of dark matter particles in the same volume that en- 
closes the stars, we would again see preferential loss of dark 
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Figure 3. In the first, second, and third row we plot the bound fraction, measured at 2.5 Gyr of harassment, of dark matter (/dm); stars 
(/*), and globular clusters (/gc) respectively, against the orbit type; plunging (orange-red colours) or circular (light-dark blue colours). 
In the final row, the ratio of foc/f* ls plotted indicating if stars or globular clusters have been preferentially lost. Columns indicate the 
model galaxy, as indicated by the column header. A particular orbit (e.g. orbit CI) has the same colour symbol in every panel. Circular 
symbols are results from the harassment simulations. A star symbol indicates a no-harasser orbit, where mass loss is purely from the 
tidal forces of the smooth background potential of the cluster. Please see Section [32] for a detailed discussion. 



matter (as shown in lSmith et alj|2010l ). This indicates that 
the orbits of the dark matter particles are more plunging and 
eccentric than the stars - even if the halo is tidally truncated 
far beyond the stars, dark matter amongst the stars can still 
be effected. Our circular orbits (orbits C1-C8) tend to lose 
more dark matter than our plunging orbits (orbits P1-P8). 
However the stochastic nature of harassment leads to some 
(especially violent) plunging orbits losing more dark matter 
than some (especially gentle) circular orbits. For any specific 
orbit, the amount of dark matter is fairly insensitive to the 
galaxy model. This can be seen by the very similar distri- 
bution of points along the top row. For example, increasing 
halo concentration from c=5 to c=30 can reduce dark mat- 
ter losses by a modest ~ 15%. The 'small disk' model suffers 
the least dark matter losses, due to the tight gravitational 
binding between stars and local dark matter. 

Now comparing top row to second and third row, it can 
be seen that baryonic losses only become significant once 
/dm ~ 0.15 (as indicated by a horizontal dashed line in the 
top row). We shall better quantify the threshold dark mat- 
ter fraction when baryons begin to be stripped in Section 
13.41 However as a result, our plunging orbits typically do 
not result in significant baryonic mass losses - hence the 
clustering of points about high bound fractions. Two excep- 
tions are the same two (especially violent) plunging orbits 



seen in the top row, where baryonic losses are roughly equal 
to the circular orbits. There is clearly large scatter in the 
baryonic mass loss for circular orbits. This indicates that 
once /dm < 0.15 and the baryons begin to be stripped, the 
amount of baryons that are lost is very stochastic - they 
are not sensitive to galaxy model, and depend more on the 
detailed interaction history. The one exception is the small 
disk model which is generally more robust to the effects of 
harassment in both dark matter and baryonic losses. 

The fourth row indicates that in almost every case GCs 
and stars are lost in roughly equal amounts for most mod- 
els (/gc//* ~ 1) when harassment is weak (for most of 
the plunging orbits). However when harassment is stronger, 
GCs are preferentially lost. This is related to their relative 
concentration - half the mass of stars can be found within 
r~1.3 kpc, whereas half the mass of GCs are found within 
r~2.9 kpc. When the GCS is extended beyond the stars 
(the 'ext GCS' model), the GCs are preferentially lost to 
an even greater extent. Only when the GCS is very con- 
centrated (the 'cone GCS' model), are stars preferentially 
lost. This is a natural e xpectation in the framework of the 
impulse approximation (ISpitzerl 1 19581 ; iGnedin et all Il999l ; 
iGonzalez-Garcfa et al.ll2005l ). 

We now consider the 'no-harasser' simulations, where 
harasser galaxies have been removed from the cluster model. 
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Except for the removal of the harasser galaxies, these are the 
same plunging and circular orbit as previously used in the 
harassment simulation. However, now mass loss is purely 
from the tidal forces of the smooth background potential 
of the cluster. It can be clearly seen that the star symbol 
(indicating a no-harasser simulation) appears above all the 
circular symbols (indicating the harassment simulations) in 
the upper, second and third row of Figure [3] This indicates 
that in all cases, the total mass loss from the smooth poten- 
tial alone sets a lower limit on the mass loss when harassers 
are included. Therefore the addition of the harasser galaxies 
provides a source of additional mass loss beyond that of the 
background potential alone. This additional mass loss occurs 
from high speed tidal encounters, and the quantity of addi- 
tional mass lost is very stochastic as it depends on the num- 
ber of close encounters along that particular orbit. This is 
valid to first order, but we add that the inclusion of harasser 
galaxies does also modify the orbit of the model galaxy. In 
this way, harassers also cause second order changes in the 
mass loss from the background potential as well. 

In summary, we have covered a considerable range of 
luminous early-type dwarf properties. These models should 
encompass the properties of the majority of the luminous 
early type dwarf galaxies, which contain sufficient numbers 
of GCs to be useful for enclosed mass measurements. We 
subject them to a wide range of strength of harassment. We 
find that their mass loss is rather insensitive to the param- 
eters we have considered. Baryonic mass loss begins only 
when the bound dark matter fraction is low for all models. 
Once baryonic mass loss begins, the amount that is lost is 
very stochastic, and does not seem to depend sensitively on 
the galaxy parameters. The one exception is the model with 
the small disk, who is more robust to harassment in general. 
For all models, GCs are tidally stripped preferentially over 
stars, with the single exception of the model with a very 
concentrated spatial distribution of GCs. 

3.3 Specific frequencies: Comparison with 
observations 

In Figure [4] we present a direct comparison between obser- 
vations and three of our models at the end of the simulation. 
Grey circles show the relation between the metal-poor GC 
mass specific frequencies (Tjv.wue) of actual Virgo dEs (from 
Peng et al. 2008 and Miller & Lotz 2007) and the host galaxy 
stellar mass. Our initial model has M* = 3 x 10 9 Mq and 
60 GCs - all assumed to be m etal-poor for the sake of this 
comparison (jPeng et al.ll2006T l - and therefore Tjv,m u <3 = 20. 
Different models and orbits are coded with different sym- 
bols and colours, respectively. Thus, plunging orbits are rep- 
resented in red, while blue symbols correspond to circular 
orbits. Stars, squares and triangles indicate the final stages 
of the standard, the 'extended GCS' and the 'concentrated 
GCS' models, respectively. It is obvious that both the inter- 
nal configuration of galaxy components and its orbit within 
the cluster play a relevant role. Our circular orbits result 
in the most significant mass losses, both for st ars and GCs. 
Our p lunging orbits inflict less damage (also see lSmith et alJ 
|2010T ). and it is tempting to suggest they are more probable 
for actual dEs harbouring rich GCSs (we shall study this 
further in a follow-up investigation). However, the initial 
GCS spatial distribution also influences the final location of 
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Figure 4. The metal-poor GC mass specific frequency (Tv blue.) 
versus stellar mass relation of actual (grey circles) and simulated 
Virgo dEs. The initial model consists ofaM t = 3x 10 10 M Q dE 
with Tjv blue — 20. Different galaxy models and orbits are coded 
with different symbols and colours, respectively. For any given 
galaxy model, circular orbits (blue) preferentially result in much 
larger mass losses than the corresponding plunging orbits (red). 
The initial spatial distribution of the GCS also plays a relevant 
role, as expected from the impulse approximation. Thus, models 
with more extended GCSs (squares) tend to preferentially lose 
GCs with respect to stars - therefore decreasing their original 
Tjv,6iue- The opposite is true for models with more concentrated 
GCSs (triangles), while the standard models (stars) lose both 
components with a similar efficiency. 

the remnant in this diagram. Models having more extended 
GCSs (squares) tend to lose a larger fraction of GCs than 
stars, and therefore their TM,biue systematically decreases 
from its initial value (c.f. Figure [3]). The standard model 
(stars) shows an almost flat behaviour, indicating that both 
GCs and stars are lost with a similar efficiency. Finally, only 
models with concentrated GCSs (triangles) exhibit a Tjv,6iue 
increase, and this is always relatively moderate - a recent, 
harassment-induced origin from more massive progenitors 
appears to be unfeasible for the lowest-mass, highest-Tjv,M ue 
Virgo dEs. These results are consistent with simple theoret- 
ical considerations of interactions in the impulsive regime 
(see ISanchez-Janssen fc AguerrJlini for a detailed discus- 
sion) . 

3.4 The amount of dark matter remaining when 
baryons begin to be stripped 

We re-emphasise that the stars begin to be stripped when 
the bound dark matter fraction /dm becomes low. To better 
quantify the required value of /dm, we record the value of 
/dm when /* reaches 0.9. In other words, when 10% of the 
stars have been stripped, we note the bound dark matter 
fraction. This is denoted /dm(/* = 0.9), and is plotted on 
the y-axis of Figure [5] (black crosses). Along the x-axis is the 
galaxy model. For every model galaxy, we obtain the average 
and standard deviation of the value of /dm(/* = 0.9) for all 
16 orbits. Despite a wide range in strength of harassment, 
and galaxy model properties, the /dm(/* = 0.9) of bound 
dark matter particles (cross symbols; black) is roughly con- 
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Figure 5. Fraction of bound dark matter particles remaining 
when the bound fraction of stars reaches 0.9 (when 10% of the 
stars have been unbound) for each galaxy model. Solid line & cross 
symbols (black) is the bound dark matter fraction. Dashed lines 
are fractions of the pre-harassment density of the dark matter 
measured within a spherical volume of radius 15 kpc (circle sym- 
bols; red) or radius 7.5 kpc (square symbols; green). Significant 
loss of stars occurs when the bound dark matter fraction reaches 
about 10 — 15% and this value is fairly model independent. This 
is when the dark matter density surrounding the globular clusters 
has roughly halved. 



stant with a value between 0.15-0.20. In other words, once 
the remaining bound dark matter fraction has fallen below 
0.15-0.20, significant stellar stripping begins. This result can 
be seen, by eye, for the two orbits shown in Figure [2l but 
now we see it remains valid for a range of galaxy models 
and over a much larger variety of orbits. Similar results are 
seen in simulations of dwarf galaxies under going mass loss as 
they orbit within the Milky Way potential |Penarrubia et all 
120081 ) , despite the fact that mass loss is often driven by high 
speed galaxy-galaxy encounters in our models. 

We also measure the total number of dark matter parti- 
cles within spheres of radius 15 kpc and 7.5 kpc, centred on 
the galaxy centre, when 10% of the stars have been stripped. 
We normalise these values by their initial value, prior to ha- 
rassment, and plot these in Figure [S] as well (red circles, 
green squares). When 10% of the stars had been stripped, 
typically the dark matter density surrounding the stellar 
disk and GCs had fallen by roughly a factor of 2, although 
the value is quite scattered between models. This explains 
why stars and GCs remain mildly affected, even when ~ 80% 
of the total dark matter has been unbound - the density of 
dark matter surrounding the stars and GCs has only halved. 

In summary, stars and GCs begin to be stripped in sig- 
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Time (Gyr) 

Figure 6. Upper panel: The evolution of the 20, 40, 60, 80, 95% 
Lagrangian radii of the GCs, for the standard model on the orbit 
P4 trajectory. While the inner most GCs are unaffected, some 
GCs are stripped, while some take up a more extended configu- 
ration. This is mainly in response to the reduced potential well in 
which the GCs orbit as mass is stripped from the galaxy. Lower 
panel: The evolution of the virial ratio of the bound GCs shows no 
significant deviation from virialised when outer GCs are stripped, 
and demonstrates the new more extended configuration of GCs 
quickly returns to a dynamically relaxed state within the new 
potential well surrounding the GCs. 



nificant quantities (> 10% stripped) once the total bound 
dark matter fraction /dm falls below 10-15%. This is true for 
all models, and so is insensitive to the the fairly large range 
of galaxy parameters we have considered. In the following 
section, we investigate to what extent the harassment and 
stripping of GCs impacts their dynamical state. 



3.5 The dynamical state of the globular cluster 
system (GCS) 

3. 5. 1 Evolution of the spatial distribution of the GCS 

The GCS of late-type d warfs typically do no t extend be- 
yond ~2 effective radii (|Georgiev et alj l2008h . Yet in dEs 
the GCSs appear t o extend to 4-7 gala xy effective radii 
(|Beaslev et all 120061 ; iBeaslev et ai"1l2009r ). Here we test if 
high speed tidal encounters could result in a more extended 
spatial distribution of the GCs that remain bound to a dwarf 
galaxy. 

In Figure[S]we plot the time evolution of the Lagrangian 
radii, and virial ratio of the GCS of the standard model, 
for orbit P4. This trajectory results in fairly mild harassment 
with mild star and GC losses (< 10%), but still substantial 
dark matter losses (~ 75%). The upper panel shows the 20, 
40, 60, 80, & 95% Lagrangian radii of the GCs (from bottom 
to top respectively). The innermost GCs are unaffected even 
when the outer GCs are stripped. However some of the GCs, 
that are still bound, take up a new, more extended distribu- 
tion that appears to be stable with time. This is primarily 
in response to the reduced potential well in which the GCs 



2 The virial ratio is defined as the absolute value of the total 
kinetic energy of the GCs divided by the total potential energy 
of the GCs. Thus a virialised system has Q=0.5. 
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Figure 7. The fractional amount that a radius containing half 
the bound GCs increases by the end of the simulation. This is an 
indicator of how the bound GCS has expanded spatially, and is 
plotted against the bound dark matter fraction. Symbols indicate 
the number of bound GCs after harassment (as shown in the key). 
There is no clear difference in behaviour between galaxy models. 
Typically the spatial distribution of their GCSs experience no 
more than a factor of two increase, due to harassment. However 
as the galaxies approach complete disruption (/dm < 5%), the 
remaining bound GCs (Nqq < 25) may expand by as much as a 
factor of 4-8. 



orbit. The potential well is reduced by harassment induced 
mass loss. The lower panel shows the evolving virial ratio of 
the bound GCs. This reveals that the new extended GCS 
is indeed stable, and has quickly relaxed back into dynami- 
cal equilibrium with the more shallow potential well within 
which the GCs orbit. 

To test the extent to which this 'puffing-up' of the 
bound GCS can occur, we measure the fraction by which 
the 50% Lagrangian radius of the bound GCs changes for 
the standard model after 2.5 Gyr of evolution, for all 16 or- 
bits. For example, if the fraction is 2, then the bound GCs 
have roughly doubled their extension. Of course unbound 
GCs can be stripped and moved to much larger radii. How- 
ever we emphasise that here we consider the change in the 
spatial distribution of the bound GCs, as these will remain 
in their new extended distribution. We plot this against 
the final bound dark matter fraction /dm of all models, 
along all orbits, in Figure [7] We find no clear trend with 
model galaxy, therefore symbols represent the final remain- 
ing number of bound GCs (as illustrated in the key). Al- 
though there is some scatter, as a rule of thumb most galax- 
ies with /dm > 0.05 expand by as much as a factor of 2-3. 
Only when the galaxy is much closer to complete disruption 
(/dm < 0.05), can the bound GCs expand more substan- 
tially by factors as much as ~4-8. We repeat all measure- 
ments using the 75% Lagrangian radius of the bound GCs, 
instead of the 50% Lagrangian radius. We find that on av- 
erage the amount of expansion is only slightly greater (by a 
factor of 1.20±0.35) if we use the 75% Lagrangian radius. We 
defer a detailed study of the link between orbit and spatial 
distribution to a future paper. 

In summary, due to harassment the current distribution 
of bound GCs about a cluster dE might be more extended 



than the original distribution by roughly a factor of 2-3 if 
the galaxy is not close to complete disruption (/dm > 0.05). 

3.5.2 How close to virialised are harassed GCSs? 

There are many sources of uncertainty that may effect a 
measurement of the enclosed mass using GC dynamics. How- 
ever one of the most fundamental of these is the assumption 
that the GCS is in dynamical equilibrium within the poten- 
tial of the galaxy. We first consider the virial state of bound 
members of the GCSs under the influence of harassment, 
before then considering other sources of error. We empha- 
sise that we only consider the dynamics of the bound GCs 
in calculating the virial ratio. This is useful as how close the 
bound members of the GCS are to dynamical equilibrium 
sets a lower limit on the uncertainty in any enclosed mass 
estimate, with other sources of error only adding additional 
uncertainty. 

In Figure [5] the virial state of the bound GCs of each 
model, after 2.5 Gyr of harassment, is plotted against the 
final bound dark matter fraction /dm- Symbols represent 
the final number of bound GCs of that particular model 
(see key for reference). We include all model galaxies, for 
all 16 orbits with the exception of those who have no GCs 
remaining with which to measure the virial ratio. The trend 
has no clear dependency on galaxy model. 

On average the virial ratio is 0.53±0.11 for all models 
with /dm >0.03. The mean (dashed horizontal line) and 
standard deviation (dotted horizontal line) are shown on 
the plot. If /dm >0.03, the full range of virial ratio for our 
models is Q=0.4-0.7. For the upper-bound value of Q=0.7, 
the enclosed mass would be overestimated by at least 40% 
if the GCS were assumed to be in dynamical equilibrium. 
To put it another way, mass estimates are - at most - 40% 
overestimated by assuming virial equilibrium, even when a 
galaxy has lost 97% of its dark matter through harassment, 
and dark matter density surrounding the GCs has fallen by 
more than a factor of two! 

However, once /dm falls below the critical value of 
/dm=0.03, the spread in virial ratio of the GCS grows 
(Q=0.19-0.99). If virial equilibrium is assumed in this 
regime, then enclosed mass estimates could be over- or 
under-estimated by a factor of ~ 2 in the most extreme 
cases. Some of this increased scatter can be attributed to 
the low-number statistics of measuring the virial ratio with 
< 10 bound GCs. However we note that there is assymetry 
in the distribution of Q values for /dm <0.03, with 57% 
found at Q>0.6 compared to only 11% found at Q<0.4. This 
would not result from simple low-number statistics, and in- 
deed some of the Q>0.6 points have between 15-19 GCs. 
Therefore the tidal encounters of harassment can maintain 
the GCSs in a super-virial state (Q>0.6). 

To further confirm that a boosted virial ratio is a gen- 
uine property of the GCS where the galaxy has /dm < 0.03, 
and not just due to the low number of GCs remaining, we re- 
peat three of the standard model harassment orbits, except 
we replace the 60 GCs of the standard model with 6000 
GCs (but with the same total mass in GCs). We refer to 
this model as the 'high N' model. We choose three orbits 
that result in /dm < 0.03; orbit CI, C5 and C8 (see Table 
[T]for remaining bound fractions with the standard model). 
Clearly this number of GCs is highly unrealistic for a galaxy 
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Figure 8. The virial state of the bound GCs of the GCS, at the 
end of the simulation as a function of the bound dark matter frac- 
tion. Symbols indicate the number of GCs that remain bound (in- 
dicated in the key) . We see no obvious dependency on the galaxy 
model used. When the bound dark matter fraction is very low 
(/dm < 0.03) and the galaxy is close to complete disruption, the 
virial state of the bound GCs may become super-virial (Q>0.6). 
In this regime, the assumption of virial equilibrium could result in 
errors of as much as a factor of 2. However, when /dm > 0.03, the 
virial state is typically close to virialised: Q=0.53±0.11 (dashed 
line is mean, dotted lines are standard deviation). The upper- 
limit is Q=0.7. Therefore, if a virialised GCS is assumed, the 
total enclosed mass will be under- or over-estimated by (at most) 
40%. 



of this luminosity, but the virial ratio of the 'high N' model is 
free from low-number effects. Bound dark matter and stellar 
fractions are identical. As we have only changed the number 
of GCs, this would be expected. Bound GC fractions differ 
by no more than 10%, and therefore are not sensitive to the 
number of GCs we initially assume. We calculate the aver- 
age virial ratio of the standard and high N model over the 
2.5 Gyr of harassment. Orbit CI has Q=0.58±0.11 (stan- 
dard model) and Q=0.60±0.07 (high N model). Orbit C5 has 
Q=0.65±0.08 (standard model) and Q=0.65±0.06 (high N 
model). Orbit C8 has Q=0.58±0.09 (standard model) and 
Q=0.59±0.11 (high N model). We see that the average Q 
is always boosted to be supervirial, independently of GC 
number, but the standard deviation of the standard model 
Q is greater due to low number effects. This underlines that 
when /dm < 0.03, dynamically the GCSs are maintained 
supervirial by tidal encounters. 

In summary, galaxies that have lost up to 97% of their 
dark matter can harbour a GCS whose virial state may cause 
enclosed mass estimates to be under- or over-estimated by no 
more than 40%. Therefore deviations from virial equilibrium 
play only a minor role in the mass estimate uncertainties. 
However once < 3% of the bound dark matter remains, the 
bound members of the GCS may become super-virial as the 
galaxy approaches complete disruption. In this < 3% regime, 
the assumption of virial equilibrium could result in mass 
estimates being incorrect by as much as a factor of 2. It 
is difficult to see how the virial state of the GCS could be 
determined independently. There are also additional sources 
of error beyond the assumption of virial equilibrium, that we 
shall consider in the following section. Therefore these errors 



can be considered lower limits on the true uncertainties in 
an enclosed mass estimate. 



3.5.3 Additional errors affecting pointed observations 

To convert the dynamics of a GCS into an enclosed mass, 
and in the absence of rotation, an equation may be used of 
the following form: 



M obs (< r i GCS ) = n r 



t GCS 



(6) 



where M b s (< ^Igcs) ls the total mass enclosed within a 
projected radius containing half the clusters of the GCS 
(r i GCS ), n is a normalising constant that differs depending 
on the mass distribution, and a is the line-of-sight velocity 
dispersion of the GCs within ri GCS . 

In Equation [6l the GCS is assumed to be in dynamical 
equilibrium. As demonstrated in Figure|8l the GCS may not 
be in dynamical equilibrium, causing a first source of error. 
Also the GC and mass distribution is assumed to be spher- 
ical, and the observed velocity dispersion may be assumed 
to be isotropic. In a pointed observation of the GCS of a 
galaxy orbiting within the cluster environment, a number of 
these assumptions may not be valid. We refer to these issues 
as 'line-of-sight' effects - additional sources of error due to 
the fact that the galaxy is observed in projection. For ex- 
ample, the GC distribution may not be exactly spherical at 
any instant - even in isolation - due to their low numbers. 
Furthermore, along its orbit a galaxy may be stretched at 
pericentre and contracted at apocentre, by the tidal field of 
the cluster. This may further exasperate how non-spherical 
the GC distribution is, as well as how isotropic their ve- 
locity dispersion is. Finally, along the line-of-sight, high ve- 
locity unbound GCs may appear close to the body of the 
galaxy, and be included in the sample, potentially resulting 
in boosting of the velocity dispersion (and an overestimate 
of the enclosed mass). 

To use Equation [6] we must first find rj which is den- 
sity distribution dependent. For a simple mass distribution, 
7] can be calculated analytically. The mass distribution en- 
closed by the GCS is a combination of the inner NFW dark 
matter halo, and an exponential disk of stars. For simplic- 
ity, we numerically calculate rj for each of our galaxy models: 
While evolving our model galaxies in isolation, we measure 
r i GCS , the velocity dispersion of the GCS within ri GCS , and 
finally the total mass (in dark matter and stars) enclosed in 
ri GGS . We can then solve Equation [6] for rj. For the standard 

model we find r\ = (5.8 ± 1.7) x 10 5 M kpc~ x (km s -1 ) -2 . 
We also calculate 77 for the other galaxy models - it does not 
vary substantially between the models and the average value 
is r, = (5.7±1.7) x 10 5 M kpc" 1 (km s" 1 )" 2 (±30% errors). 
Therefore mean enclosed masses calculated, using Equation 
have 30% errors from n alone, as well as additional errors 
from 'line-of -sight' effects and the assumption of dynamical 
equilibrium. 

To quantify how much these combined effects may in- 
fluence our mass estimates, we first use Equation [6] on 
all our model galaxies after 2.5 Gyr of harassment. For 
this, we use the average value found for rj (rj = 5.7 x 10 5 
M@ kpc~ x (km s -1 ) -2 ) and, for each individual model, we 
choose three independent lines of sight (along the x, y, and z 
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Figure 9. The observed mass over the real mass as a function 
of the bound dark matter fraction. The observed mass M t, s 
is calculated using Equation [6] whereas M rea i is the real mass 
measured directly from the simulation. If /dm is greater than 
3%, the observed mass can be estimated to within a factor of 
two due to a combination of deviation from dynamical equilib- 
rium, uncertainty in r\ (see Equation [6}, and additional line- 
of-sight effects (such as unbound GCs). The dashed and dot- 
ted horizontal lines show the average and standard deviation 
(M obs /M rea i=1.05±0.41) for / DM > 3%. Symbols indicate the 
number of GCs found within ri GCS , and therefore used to cal- 
culate the observed mass. In some cases, the mass is measured 
to within a factor of two using only < 15 GCs. However for 
/dm < 3%, when the galaxy is close to complete disruption, mass 
estimates can be substantially overestimated. 



axis). This gives us 384 data points in total (3 lines of sight 
and 128 models). We exclude galaxies that are completely 
disrupted (/dm = /*= 0.00) leaving 330 data points in total. 
In order to mimic the observations, we assume a field-of-view 
that includes GCs up to 8 effective radii (10 kpc) from the 
galaxy. We later demonstrate our results are not sensitive 
to this rather arbitary choice of field-of-view. We then mea- 
sure ri GCS of the GCSs in the field-of-view, and the velocity 
dispersion of the GCs within r i GCS . Then we calculate the 
estimated enclosed mass for each line-of-sight (M b a ) using 
Equation [5] Finally we divide each value by the true en- 
closed mass within a radius equal to ri GCS along that line 
of sight (M rea i), providing us with the quantity M bs/M rea i. 
This ratio directly tells us by what factor we measure the 
true enclosed mass, and is plotted on the y-axis of Figure [5] 
We do not differentiate between different galaxy models in 
Figure[9] as we find no obvious dependency on galaxy model. 
Instead we provide each data point with a symbol according 
to the number of GCs found within ri GCS (and thus the 
number of GCs used to calculate the velocity dispersion for 
that data point). 

Below the critical bound dark matter fraction /dm = 
0.03, the mass estimates typically become overestimated. Of 
all galaxies in this regime, 35% have their real mass over es- 
timated by more than a factor of two. Only 17% recover the 
true mass to within 20%. As shown in Figure [8] the bound 
GCs are often supervirial when /dm > 0.03. However this 
could only result in mass over estimates by a factor of two of 
the true mass at most. The key cause of such boosted mass 



estimates appears to be primarily unbound GCs, seen down 
the line of sight. With such a low dark matter fraction, the 
galaxy is close to complete disruption, and a large number of 
the previously bound GCs may suddenly find themselves un- 
bound. Therefore a large number of unbound GCs are seen 
in projection to be close to the galaxy, and their velocity 
dispersion is large, resulting in boosted mass estimates. 

For /dm > 0.03, the observed mass recovers the true 
mass well for all data points, but with significant scatter - 
M b B /M Taa i is on average 1.05T0.41 (the horizontal dashed 
line is the mean, and the horizontal dotted lines are the 
standard deviation). The scatter is a little larger than the 
±20% standard deviation we get from the virial ratio alone, 
due to the additional line-of-sight errors and uncertainty in 
the normalising constant r\ (see Equation [6]) . However, tak- 
ing into account the full range in M ^ B /M Tea i seen, above 
/dm ~ 3%, mass estimates are accurate to within a fac- 
tor of 2. We note that some green diamond symbols can be 
found close to M b s /A/ roa i = 1. Therefore, mass measure- 
ments that are accurate to within a factor of two can be 
made using the velocities of only 10-14 GCs. 

We repeat all measurements using a field-of-view en- 
compassing all GCs within 12 effective radii (15 kpc) of each 
model, and find our statistics are in close agreement to when 
our field-of-view was 8 effective radii (10 kpc). Thus the ac- 
curacy of our mass estimates are not sensitive to our rather 
arbitary choice of field-of-view. Unbound GCs do not sig- 
nificantly increase enclosed mass errors in the /dm > 0.03 
regime as once they are unbound, they quickly separate from 
the main body of the dwarf galaxy and few are seen along the 
line-of-sight. This remains the case even when we increase 
our field-of-view out to 12 effective radii. These stripped 
GCs can make an important contribution to the intraclus- 
ter component, as seen recently in nearby clusters l|Lee et alj 
l20irj IPeng et al.ll201ll ; IWest et alj|201lh . We defer a study 
of the properties of the stripped GCs for a future paper. 

In summary, as long as the remaining bound dark mat- 
ter fraction /dm > 0.03, and assuming a non-rotating GCS, 
an enclosed mass may be calculated using Equation [6] and 
77 = 5.7 x 10 5 M kpc^ 1 (km s -1 ) -2 . This method of mass- 
measurement can be applied to galaxies with as few as 10-15 
GCs. However non-virialised GCSs, line-of-sight effects, and 
uncertainty in the assumed rj provide sources of error. As a 
result, we find the calculated enclosed mass agrees with the 
true enclosed mass to within a factor of two if the bound 
dark matter fraction is above a critical value /dm > 0.03. 



3.6 Determining if a galaxy is above the critical 
dark matter fraction 

Using GC dynamics to measure the total enclosed mass is 
clearly only reliable if a galaxy has more than the critical 
amount of bound dark matter (/dm > 0.03). One possible 
means to test if this is the case is by measuring the dynamics 
of the stellar disk. 

It is well established that, if harassment is strong, it 
can effectively heat fragile stellar disks and convert them 
from rotationally supported systems into dispersion sup- 
ported systems. Therefore we test the hypothesis that the 
amount of rotation might provide constraints on the bound 
dark matter fraction. 

For all our model galaxies, falling along all orbits, we 



Harassment of globular cluster systems 13 




+ Mastropietro (2005) 
9 These simulations 
I i I 



0.5 



bound dark matter fraction f, 



DM 



Figure 10. The ratio of the final to initial i> pca k (the peak in 
the rotation curve) as a function of the dark matter fraction. 
Filled circular symbols (black) with error bars are the model 
galaxies. Plus symbols ( red) are harassed dwarf galaxies from 
iMastropietro et al.l i 20051) . and represent lower limits. The peak 
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as dark matter is lost. As /dm falls below 0.05, the peak rotation 
velocity falls rapidly. ^ pea k i s typi can y reduced to 0-20% of its 
prc-harassment value as /dm approaches zero. 



produce the rotation curve of the stellar disk after 2.5 Gyr 
of harassment, and record the peak value of rotation %, ea k. 
To do so, we choose our lines of sight to be along the plane 
of the disk, so as we do not have to correct for inclination. 
We then calculate the average line of sight velocity through 
the disk, in bins of distance from the disk centre. Our bins 
are 1 kpc wide, and we typically measure the rotation curve 
out to 10 kpc which is sufficient to encapsulate u pe ak- We 
reject bins with < 50 particles to reduce noise in rotation 
curves at large distances from the disk centre. Although our 
galaxy disks are dispersion dominated from the beginning, 
this method clearly and reliably recovers the rotation curve. 

We note the w pea k with estimated errors for each or- 
bit, and divide these by the initial « pea k before harassment. 
This gives us the fraction of the initial v pca k as a result 
of harassment. We plot this fraction against the remain- 
ing bound dark matter fraction /dm in Figure [TO] Our 
simulations are shown as filled circles (black) with error 
bars. We also include da ta points from the simulations of 
IMastropietro et al.l (|2005l ) (red plus symbols), where we take 
v P cnk values from their provided rotation curves. The ini- 
tial u pea k of their model galaxy is not specifically given, 
so we assume f pea k = 92 km s , based on their Fig. 1 in 
IMastropietro et al.l (|2005l ). Also the inclination of their disks 
is unknown, therefore all of their data points can be consid- 
ered lower limits. 

The peak value in the rotation curve of harassed dwarfs, 
Vpc&k, tends to decrease as dark matter is lost. At a fixed 
/dm, there is clearly some variation in the fraction of the 
initial u pca k- However the upper-most value of the fraction 
falls rapidly when /dm falls below about 0.1. w pca k is typi- 



cally reduced to 0-20% of its pre-harassment value as /dm 
approaches zero. In short, the amount of rotation in the stel- 
lar disk is strongly reduced when the remaining bound dark 
matter becomes small. 

We therefore suggest that rotation in the stellar disk 
could provide an indication of whether GCS dynamics are 
useful to estimate the enclosed mass. For example, consider a 
disk that is strongly supported by rotation. It is likely that 
such a galaxy still contains sufficient bound dark matter, 
that its GCS dynamics are useful to estimate the enclosed 
mass. Otherwise the disk would have had to be considerably 
more dominated by rotation, prior to mass loss, than it is 
now. However, we emphasise that lack of rotation does not 
necessarily imply low dark matter content - the disk may 
well have been dispersion dominated before infall into the 
cluster. 



3.7 Can the mass loss by harassment be detected? 

A clear finding in all harassment simulations is that mass 
loss is substantial and unavoidable. Even the 'weakly' ha- 
rassed galaxy models of our sample lose over half their dark 
matter. This poses a question - is this mass loss detectable 
through observations of GCS dynamics? Perhaps by com- 
paring a sufficient sample of field and cluster galaxies, and 
looking for an offset in enclosed mass? 

We attempt to answer this question by calculating the 
ratio of the 'observed' mass (calculated using Equation [6| 
over the measured, pre-harassment mass enclosed in the 
same ri GCS of the GCS. We plot this against the bound 
dark matter fraction in Figure [TT] along three independent 
lines of sight, for all of our model galaxies along all 16 or- 
bits (384 data points in total). We exclude galaxies that are 
completely disrupted (/dm = /*= 0.00) leaving 330 data 
points in total. The non-virialised state of the GCS, com- 
bined with 'line-of-sight' errors, can result in a fraction of 
the galaxies appearing to contain equal quantities (or more) 
of dark matter before and after harassment. We note that it 
is not possible for the galaxies to truly have more dark mat- 
ter than before infall - they only appear to have increased 
mass due to their GCs being supervirial, and due to un- 
bound GCs and other projection effects. In the following we 
shall split our sample into galaxies that appear unaffected 
by harassment, and those which appear to have been signifi- 
cantly affected. We emphasise the word 'appear' to be clear 
that this is observed mass we would measure using Equation 
[6] and not necessarily their true enclosed mass (e.g. see Fig- 
ure [9]). The dotted lines on Figure [TT1 indicate galaxies such 
as these, who appear almost unaffected by harassment, with 
an estimated enclosed mass of 0.8-1.2 of their initial mass. 
All galaxies below the dashed line appear heavily affected, 
as they appear to have lost more than half their initial mass. 

The relative number of galaxies which appear affected 
or appear unaffected is a function of the bound dark mat- 
ter fraction that the galaxies in the sample have. We note 
that the GCS dynamics only probes the enclosed mass sur- 
rounding the GCs. Even when 85% of the dark matter halo 
has been stripped (/dm=0.15), typically the mass of dark 
matter surrounding the GCs has only fallen by roughly a 
factor of 2 (see Figure [5]), where as the mass loss of stars is 
almost negligible (/*=0.9). Therefore measurements of the 
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Figure 11. The ratio of the 'observed' mass over the measured 
pre-harassment mass (measured within i"Xqqq)^ plotted against 
the bound dark matter fraction. Some galaxies appear to contain 
equal dark matter mass as prior to harassment (M b B /Mi n i t =l). 
The dotted horizontal lines arc placed at M o b s /Mi n i t =0.8 and 1.2, 
while the dashed horizontal line is at M o b s /M; n ; t =0.5. For /dm > 
0.03, where enclosed masses are well measured, the fraction of 
galaxies that appear unaffected by harassment (i.e. between the 
dotted lines) is roughly equal to galaxies that show more than 
half their mass has been lost (i.e. below the dashed line). 



enclosed mass would not be expected to show a substantial 
mass deficit until /dm < 0.15. 

If we consider only galaxies with /dm < 0.03, galax- 
ies that appear affected by harassment (74%) outnumber 
galaxies that appear unaffected (4%) considerably. However 
if we consider only galaxies with /dm > 0.03, where the en- 
closed mass is well measured, galaxies that appear affected 
by harassment (23%) are found in roughly equal numbers to 
galaxies that appear unaffected (26%). 

In summary, in a sample of cluster dwarf elliptical galax- 
ies, the number of galaxies that appear unaffected by ha- 
rassment, compared to the number that appear to have lost 
half their mass, depends on the bound dark matter frac- 
tions of the galaxies in the sample. In a real galaxy sample, 
the bound dark matter fraction of individual galaxies might 
be expected to vary considerably, depending primarily on 
orbit but with significant scatter due to the stochastic oc- 
currence of high speed galaxy-galaxy encounters. Galaxies 
that appear affected significantly outnumber galaxies that 
appear unaffected if there is little dark matter remaining, 
and galaxies are close to complete disruption. However if the 
galaxies have > 3% of their dark matter remaining, galaxies 
that appear affected are found in roughly equal numbers to 
galaxies that appear unaffected. 



4 SUMMARY & CONCLUSIONS 

We conduct N-body simulations of cluster early-type dwarf 
galaxies, surrounded by a globular cluster system (GCS). 
Our standard model is an approximate model for luminous 
Virgo early- type dwarfs, and we consider galaxy models with 
a range of halo mass, halo concentration, disk size, and spa- 
tial distribution of the GCS. We expect that the range of pa- 
rameter space we have covered should encompass the prop- 



erties of the majority of the early type dwarfs, who contain 
sufficient numbers of GCs to be useful for mass measure- 
ments. We evolve each galaxy model for 2.5 Gyr within a dy- 
namic, time-evolving potential field that mimics the effects 
of harassment in a galaxy cluster like Virgo. Each galaxy 
model experiences widely ranging strengths of harassment. 

We first study the amount of dark matter, stars, and 
globular clusters that are stripped by harassment. We find 
that the stars and globular clusters are not stripped in sig- 
nificant quantities until the remaining bound dark matter 
fraction falls to 10-15%, and when the dark matter density 
surrounding the disk has been reduced by roughly a factor 
of two. 

We also study the effects of harassment on the dynam- 
ics of the GCS. We find that the bound GCs are typically 
close to virialised (Q=0.53±0.11) within the galaxy poten- 
tial while the bound dark matter fraction is above 3%. We 
also provide a normalising constant (77), valid when there 
is no rotational support, that enables the enclosed mass of 
dwarf galaxies to be calculated from their observed GCS 
velocity dispersion. 77 has intrinsic uncertainty of 30%. In- 
cluding additional sources of error, that will plague pointed 
observations of dwarf galaxy GCSs (e.g. unbound GCs, as- 
symetric spatial distribution of GCs, and other line-of-sight 
effects), we find that enclosed masses may be accurate to 
within a factor of 2 of their true masses. This is valid while 
the bound dark matter fraction is above 3%. 

If the bound dark matter fraction is below this critical 
3% value, the galaxy is close to complete disruption. At this 
stage, bound members of the GCS can become super-virial. 
Unbound GCs can also be found close to the body of the 
galaxy resulting in a boosted velocity dispersion and there- 
fore the GCS dynamics cannot be used as a reliable measure 
of the enclosed mass. It is therefore critical that we can de- 
termine if a galaxy has more or less than 3% of its dark 
matter still bound. 

A possible indicator for whether the galaxy has > 3% 
dark matter bound could be measurement of rotation in the 
stellar disk. When the bound dark matter fraction is reduced 
below 3%, typically the peak in the rotation curve is reduced 
to < 20% of its pre-harassment value. Therefore a galaxy 
that currently appears to be rotationally supported is likely 
to contain sufficient dark matter, that dynamics of its GCS 
can provide an estimate of the mass enclosed, to within a 
factor of 2 of its real value. 

Our key results may be summarised as follows. 

(i) Harassment results in the stripping of dark matter, 
stars and globular clusters (GCs) from our model galaxies. 
We consider a wide range of parameters for our model galax- 
ies varying halo mass, and halo concentration, disk size, and 
globular cluster spatial distribution. We expect this param- 
eter study to encompass the properties of most early-type 
dwarfs of interest for enclosed mass measurements. We find 
that mass loss is not very sensitive to the model galaxy pa- 
rameters. 

(ii) As dark matter is stripped, the bound dark matter 
fraction /dm falls. We find that /dm is a key parameter 
controlling baryonic mass loss, GC spatial distribution and 
GC dynamics. 

(iii) Stars and GCs are not significantly stripped until 
the bound dark matter fractions is very 1ow(/dm <10-15%), 



Harassment of globular cluster systems 15 



with little dependece of parameters of the model galaxies we 
consider. 

(iv) Unbound GCs are often stripped to large radii. How- 
ever, even considering bound GCs only, the globular clus- 
ter system (GCS) can expand spatially in response to the 
changing depth of the potential in which the GCs orbit. For 
/dm > 5%, the bound GCS can expand by a factor of 2-3. 
For /dm < 5%, as the galaxy approaches total disruption, 
the bound GCS can expand by a factor of 4-8. 

(v) Despite the stripping of GCs, an observed velocity 
dispersion, along a line-of-sight to the GCS, can be used to 
estimate the enclosed mass surrounding the GCS to within 
a factor of 2. This remains valid even when /dm is as low as 
3%. 

(vi) However once /dm < 3%, the galaxy is close to to- 
tal disruption, and the observed velocity dispersion can no 
longer provide reliable mass estimates. Observing a rotation- 
ally dominated stellar disk could indicate that a galaxy has 
/dm > 3%. 

In this study, we have focussed on globular cluster sys- 
tem (GCS) dynamics as a means of measuring the enclosed 
mass. However in a future publication we shall study a num- 
ber of other related topics including harassment-induced ro- 
tation in the GCS, the link between spatial distribution and 
orbit, and the fate of stripped GCs. 

The tidal potential field of a galaxy cluster has a de- 
structive influence on cluster galaxies. Substantial loss of 
dark matter is seen in almost all harassment simulations 
( Mastropictro ct al 1 l2005l ; lAguerrT fc Gonzalez- Garci'all2009l ; 
ISmith et al.ll2010h . and is almost inevitable. One reason to 
measure the enclosed mass of cluster dwarf elliptical galaxies 
might be to try to detect this mass loss. In our models, we 
find that the mass loss may be detectable. However we add 
that, within such a sample, galaxies that appear unaffected 
by mass loss may be equally common as galaxies that appear 
to have lost half their enclosed mass. Therefore a reasonable 
number of galaxies may be required to detect the mass-loss. 

We additionally note that our findings may have im- 
plications for the assembly of larger galaxies via merging of 
smaller galaxies. In the cluster environment, dwarf galaxies 
can have ~ 85% of their dark matter stripped, but keep al- 
most all their GCs. Therefore input of dark matter through 
merging may be overestimated, if a 'pre-harassment' globu- 
lar cluster to dark matter mass ratio is assumed. 
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